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ANNALS OF MATHEMATICS. 



Vol. I. March, 1884. No. i. 

THE DETERMINATION OF THE MASS OF A PLANET FROM THE 
RELATIVE POSITION OF TWO SATELLITES. 

Bv Prof. Asaph Hall, Washington, D. C. 



One of the easiest methods of determining the mass of a planet is from the 
observations of its satellites. The two principal elements that are needed are 
the periodic time of the satellite about its primary, and the semi-major axis of 
its orbit. Thus, if Maud m are the masses of the sun and the planet, and r^ 
and X are the periodic times of the planet and the satellite, we shall have very 
nearly, 

in r^ a 8 

«Q and a being the mean distances of the planet from the sun and of the satellite 
from its primary. The above is the simplest form of this equation, and should it 
be necessary in any case to consider the elliptical figure of the planet, and its mass 
when added to that of the sun, these small corrections can be applied easily. Of 
the two principal elements the periodic time of a satellite can be found with suf- 
ficient accuracy by comparing observations of it made at distant epochs, and 
hence the error arising from this source may be considered as vanishing. But 
with respect to the semi-major axis the case is different. For the determination 
of this element we have to depend on micrometrical measurements, and several 
sources of error present themselves. Thus the value in arc of one turn of the 
screw of the micrometer must be known very exactly, since this value will fre- 
quently be multiplied by a large factor, and of course an error in the value of 
the screw will also be multiplied by the same factor. With care on the part of 
the astronomer this error may be reduced to a very small quantity ; and yet it is 
surprising how often such errors remain in micrometers that have been a long 
time in use. But there is another error in micrometrical measurements that is 
not easily got rid of, and this is the personal error of the observer. This error 
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will sometimes amount to half a second of arc, and since it is not eliminated by 
the repetition of observations, each observer will have his own value of the mean 
distance of a satellite, and hence his own mass of a planet. This personal error 
may be still greater, if, as is generally the case, we have to do with the disk of 
a planet. In this case we must bisect the disk with one of the wires, or else 
measure from the limb of the planet, and both operations are difficult and tend 
to increase the constant error of the measurement. 

On account of this difficulty Struve, the present Director of the Pulkowa 
Observatory, has proposed that observations be made of the relative positions of 
two satellites for the purpose of correcting their elements, and in this way de- 
ducing the mass of the planet. Mr. Backlund, one of the astronomers of the 
Pulkowa Observatory, has already begun a series of measurements of the rela- 
tive positions of the satellites of Jupiter with a heliometer, and it is said that 
these observations are two or three times as accurate as those made by reference 
to the limb of the planet, and moreover that the complete observation of the 

relative position of two satellites requires only a 
third part of the time that is necessary to de- 
termine the position of one of them with respect 
to the center of the planet. Mr. Backlund has 
published an outline of this method, which is 
theoretically very simple. Let P be the place 
of the planet projected on the heavens, JJ the 
pole of the heavens, 5] and ^2 the satellites, and 
denote the angles and sides of the spherical tri- 
angle, as in the figure. 
The angles of position and distances are designated generally by / and s, 
with proper subscripts, and the relative angle of position and the distance of the 
two satellites by 77 and a. Then 

angle S^PS^ =p^ ~p^, angle S^S^P=p\ — -', 

and the spherical triangle S-^PS^ gives the equations, 

sin <T sin (p\ — tt') = sin s^ sin (p^ — p^ ) 

sin a cos {p\ — tt') = sin s^ cos ^2 — cos .fj sin s,^ ^os (p^ — pi)- 

These are exact formulas, but as s-^ and s^ are always less than 12' in any real 
case that we have, we may generally put 

pl^^p^, and;r'=:;r; 
cos .fj =cos .?2 ^= i» ^"d sin 5j ^5]^, sin52=i'2. sin <T:=<T, 
and hence we get 

a sin (/i —n)=.s^ sin {p^ —p^) 

a cos (/i — ;r) = 5i — 52 cos {p^ ~Pi), 
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and from these we have 

a sin - = 5j sin/]^ — s.^ sin p.^ 

a cos ;: ^= s^ cos/j — ^2 cosp^- 
Now 5j sin/i and s-^ cos/, are the rectangular co-ordinates of the first satel- 
hte with respect to the center of the planet, and the preceding equations become 

1 = Ii - Ja- 
where x and j are the relative co-ordinates of the satellites ; a result that might 
have been anticipated, since the above assumptions reduce the work to that of 
plane co-ordinates. Should it be necessary in any case we may use the rigor- 
ous formulas; or, what is easier, we may expand them and retain terms of the 
lower orders. In the present case we have therefore 

dx = dx^ — dxc^ 

dy = dy^ -- dy^. 

The differentials on the right hand side can be expressed in terms of the differ- 
tials of the elements of the orbits of the two satellites, and of the differential co- 
efficients with respect to their elements; and since the left hand side will be 
known by comparing the computed relative position with the observed, we have 
the means of forming the equations of condition for correcting the elements of 
these orbits, and therefore for determining the mass of the planet. Thus dx-^ 
and dx^ will be of the form 

where aj,s^, Pi,e-^, N-^,J-^ are the elements of the orbit of the first satellite, and 
(221 *2> -^2' ^2> ^v Ji ^fs those of the second satellite. The equations of condi- 
tion therefore that result from this method will each contain twelve unknown 
quantities ; and it will be necessary to arrange the observations so that there may 
be no indetermination in the solution. For this purpose the satellites should be 
observed when they are on different sides of the planet as well as when they are 
on the same side, and in a variety of positions. 

In the common method of referring a satellite directly to the center of the 
planet we have in the equations of condition six unknown quantities. Now the 
labor of solving a set of equations of condition increases rapidly with the number 
of unknown quantities. If i be this number, then the number of coefficients in 
the normal equations is found by the summation of an arithmetical series to be 

2 
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The number of auxiliary quantities that have to be formed in the successive 
elimination of the unknown quantities from the normals is given in the following 
table, with their first and second differences. 

No. Ai A2 

/'2 + i — 2 



«2- 


2 


2 


P- 


2 





P- 


2 


4 


P- 


2 


10 


P- 


2 


18 



~i^o 






+ 1 


-i^i 






+ 1 


— i^2 






+ 1 


-i^l 






+ 1 


-i+A 





2 
etc. 

The law of this series is evident. The sum of i — I terms of the series is 

Sum = ^('— (^'+2) P U — i) ^ i{i—i) — 2) 
2 2 2.3 

2-3 
Hence the whole number of quantities to be formed in this solution is 

(^-1)^(^+4) ^ ^(^+3) ^ ^(^+i) (^+5) 
2.3 2 2.3 

In the case of i=--6, we shall have JJ quantities to form, and when i= 12, we 
shall have 442 quantities to form. The labor of solution is nearly six times as 
great in the new method, but this is a consideration that is not perhaps entitled 
to much weight. Again, when the number of unknown quantities is increased, 
the probable error of a single equation is generally increased, but this increase 
will be small when the number of observations is large ; and if the observations 
be well arranged, this question may be left out of consideration. On the other 
hand, the gain in accuracy seems to be substantial, and the new method prom- 
ises a real advance, especially in the case of the satellites of Jupiter and Saturn. 
Another advantage is that it is a different method of making the observations, 
and in this way it offers an independent check on the old methods, and may be 
the means of eliminating some constant errors. 



